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Abstract: A spectrum of massive graviton states is present in several recent 
theoretical models that include extra space dimensions. In some such models the 
graviton states are well separated in mass, and can be detected as resonances in 
collider experiments. The ability of the ATLAS detector at the Large Hadron 
Collider to identify such states and measure their properties is considered, in the 
case that the resonances are narrow compared to the experimental resolution. The 
discovery limits for the detection of the decay mode G — > e + e~ are derived. The 
angular distribution of the lepton pair is used to determine the spin of the interme- 
diate state. In one specific model, the resonance can be detected up to a graviton 
resonance mass of 2080 GeV, while the angular distribution favours a spin-2 hypoth- 
esis over a spin-1 hypothesis at 90% confidence for resonance masses up to 1720 GeV. 
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1. Introduction 

The recently proposed localized gravity model of Randall and Sundrum |l] has 
aroused great theoretical interest, and many possible extensions and elaborations of 
this type of theory are being discussed in the literature 0. The model was initially 
motivated because it solved the weak-Planck scale hierarchy via an exponentially sup- 
pressed warp factor in a non-factorisable geometry. Despite possible problems with 
negative tension branes and stability in some of the frameworks, these problems are 
solved in other models. For now, we ignore the details of particular models and con- 
centrate on a measurement that should apply to a broad class of the models. Such 
scenarios are characterised by the existence of a series of graviton excitations, which 
may be detectable at present and future colliders. In some cases these excitations 
may be sufficiently well spaced to be detected as individual resonances. 

This work will discuss the possibility of detecting such graviton resonances in as 
model independent a way as possible, using the scenario of as a guide. However, 
the results derived do not depend on the validity of this particular scenario, but can 
be applied to any model giving rise to narrow graviton resonances. 

In the scenario of the massive graviton excitations couple with equal strength 
to the visible sector [[|. However, the higher modes being suppressed by the falling 
parton distribution functions, we consider only the lightest mode. This does not in 
any way affect the generality of our approach, as our analysis can be applied to any 
such resonances, including the higher modes, so long as the resonances are narrow 
and sufficiently separated from the other modes. This is in contrast to studies in 
which many excitations, each with small coupling, contribute to some scattering 
process ||. 
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An event generator capable of simulating the production and decay of spin-2 res- 
onances has been developed. This generator is an extension of the HERWIG 6.1 || 
simulation package. The generated events are passed through the ATLAS fast sim- 
ulation (ATLFAST H), in order to give a realistic description of detector resolution 
and efficiency. The decay channel G — > e + e~ is chosen for study, since the ATLAS 
detector has excellent energy and angular resolution for high energy electrons, and 
the Standard Model background comes dominantly from the well-understood Drell- 
Yan process. The observable cross section for this process is considered first, and 
limits on the discovery reach of ATLAS, based simply on the ability to detect the 
resonance above the Standard Model background, are derived. The angular distri- 
bution of the decay products is then considered, in order to study the ability of the 
detector to determine the spin of the resonance. 

2. The event generator 

In the model of H], a 5-dimensional non-factorizable geometry is used, with two 
3-branes of opposite tension. A graviton Kaluza-Klein spectrum is created, with a 
scale A n = Mpie~ krc7T , where Mpi is the reduced effective 4-D Planck scale, r c is 
the compactification radius of the extra dimension, and k is a scale of the order 
of the Planck scale. The geometrical exponential factor (dubbed the 'warp factor') 
generates TeV scales from fundamental Planck scales and hence offers a solution to 
the hierarchy problem, if kr c ^12. The masses of the graviton resonances are given 
by m n = kx n e~ krc7T = x n (k/ Mp^A^ where x n are the roots of the Bessel function of 
order 1. The couplings of the massive resonances are given by 1/ A n . The properties 
of the model are determined by the ratio k/Mpi. We have chosen a value of 0.01 
(at the bottom of the range suggested in ||) for this ratio, which according to [0], 
is on the edge of 95% exclusion for a first graviton excitation mass of less than 2000 
GeV. Thus, we assign a low coupling constant to the gravitons, and hence obtain 
a conservative estimate of the production cross section. This choice leads naturally 
to narrow resonances. This test scenario is used to illustrate the potential physics 
reach. 

The implementation of the resonance is largely model independent, depending 
only on a universal coupling to the Standard Model fields. This is generally the case 
for models deriving from extra dimensions, and the effective Lagrangian is given at 
first non-trivial order in 1/A n as: 

C^—hTT^ (2.1) 

where h^ v is the spin-2 field and T^ u is the energy-momentum tensor of the Standard 
Model fields. Here, we consider only the production and measurement of the lightest 
massive graviton excitation. We note however JO, that other heavier resonances are 
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expected. In the case we consider here, they are well split in mass and hence the 
lightest resonance will be produced dominantly. 

The decay G — > e+e~ is treated in the HERWIG implementation as a 2 — > 2 
process, consisting of the two hard production subprocesses qq — > e+e~ and gg — > 
e + e~ . For the qq — > e+e~ channel we do not consider interference with the Standard 
Model Drell-Yan subprocesses as the interference is negligible for our purposes. The 
code is available from the authors on request. 

For calculating the matrix elements and the total decay width of the graviton 
resonance, we made use of the results in || || [UJ. We carried out our calculations 



both by hand and by using HELAS |TTJ incorporating subprograms with spin-2 par- 
ticles ||12|| . Our implementation of the decay widths corrects a typographical error 
in || for the decay into heavy vector bosons but this has only a small effect on the 
total decay width. 



Y(G^VV) = 5- 
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40vrA2 v " " y ' V12 3 

where, as in [6], 5 = 1/2 for identical particles and ry = m v /mQ. 

The angular distributions of the possible subprocesses, in the centre-of-mass 
frame of the resonance, are shown in Table |l|. Here 9* is the angle between the decay 
electron and the beam direction in the dilepton centre-of-mass frame. The spin-2 
distributions contrast strongly with the angular distributions resulting from spin-1 
or scalar resonances. These distributions are shown in Figure [I]. 



Process 


Distribution 


99 -> G -> e+e" 
qq — > G — > e + e~ 
qq, 99 -> e+e" 
qq, 99 -> S -> e+e- 


i - cos 4 e* 

1 - 3 cos 2 9* + 4 cos 4 9* 
1 + a cos 2 9* 
1 



Table 1: Angular distributions in graviton (G), vector (V) and scalar (S) boson production 
and decay, a = 1 in Standard Model processes. 



3. Detection of the graviton resonance 

The process pp — > G — > e + e~ must be detected above the background from the 
Standard Model processes pp — > Z/7* — » e+e- . Samples of events were generated for 
a single graviton excitation with a mass ranging from 500 GeV to 2.2 TeV. The masses 
and widths of the graviton resonances are given in Table |2|, as they are computed 
with the model parameters given above. The table also shows the cross-section times 
branching ratio a ■ B in the test model and the width of the resonance reconstructed 
after the ATLAS detector simulation. In all the cases studied, the true width of the 
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resonance is much smaller than the Gaussian experimental resolution T m and can be 
neglected. Hence the results are valid for any model in which this condition is met. 
In our test model, this is the case for k/Mpi < 0.06, at which point the width is 
approximately equal to the experimental resolution, and the cross section is 36 times 
larger than we have assumed. For larger values of k/Mpi, the very high signal cross 
sections would overwhelm the Standard Model background, making the resonance 
very easy to detect, and the spin simple to determine, even with a large width. 

For each mass point, having determined the observed resonance width T m , a 
mass window is defined as ±3r m around the graviton mass. The Drell-Yan cross- 
section inside this mass window is then calculated using the HERWIG 6.1 event 
generator 0. The number of background events expected in the mass window is 
then computed, for an integrated luminosity of 100 fb _1 , corresponding to one year 
of LHC running at design luminosity. The ATLFAST simulation of the ATLAS 
detector is used to compute the acceptance for the process, which is dominated by 
the pseudorapidity (77) coverage of the tracker, \r}\ < 2.5. Cuts are also made on the 
isolation of the electrons to account for the effect of nearby hadrons on the detection 
efficiency In addition to the cuts imposed by ATLFAST, a tracking efficiency of 
90% was applied to each electron, to allow for track reconstruction losses which are 
not included in ATLFAST. These losses are estimated in [H| to be approximately 
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m G (GeV) 


T G (GeV) 


T m (GeV) 


a-B (fb) 


500 


0.068 


3.53 


281.9 


1000 


0.141 


6.02 


11.0 


1500 


0.213 


8.13 


1.20 


1700 


0.242 


8.78 


0.57 


1800 


0.256 


9.34 


0.41 


1900 


0.270 


9.66 


0.29 


2000 


0.285 


9.80 


0.21 


2100 


0.298 


10.18 


0.15 


2200 


0.312 


10.49 


0.11 



Table 2: Masses mc and widths Tq of the simulated graviton resonances, and the Gaussian 
width of the observed resonance T m after detector effects, u • B is the cross-section times 
branching ratio in the test model. 

5% at high transverse momentum if the tracks are well isolated. Our efficiency value 
is therefore conservative. The number of Drell-Yan background events (Ng) inside 
each mass window after these cuts is given in Table ^. 



m G 


Mass window 


N s 


N B 


N min 


(a ■ B) min 


(GeV) 


(GeV) 






(fb) 




500 


±10.46 


20750 


816 


142.9 


1.941 


1000 


±18.21 


814 


65 


40.2 


0.542 


1500 


±24.37 


84.3 


11 


16.5 


0.235 


1700 


±26.53 


38.6 


5.8 


12.0 


0.178 


1800 


±27.42 


27.2 


4.3 


10.4 


0.156 


1900 


±28.29 


19.2 


3.2 


10.0 


0.152 


2000 


±28.76 


13.2 


2.3 


10.0 


0.157 


2100 


±30.55 


9.4 


1.8 


10.0 


0.159 


2200 


±31.46 


6.8 


1.4 


10.0 


0.162 



Table 3: The number of signal events in the test model, Ns, expected inside each graviton 
mass window in a run of 100 fb -1 , after detector effects; the number of Standard Model 
background events, Nb; the minimum number of signal events required to detect the reso- 
nance, iVjp m ; and the minimum production cross-section times branching ratio (a ■ B) mtn 
required to detect the resonance via e + e~ production. 

The expected signal in the test model for m,Q = 1.5 TeV can be seen in Fig- 
ure |2], showing the data expected for 100 fb -1 of integrated luminosity. The signal 
is prominent above the small Standard Model background. 
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Figure 2: The number of events per 4 GeV mass bin from a graviton resonance, with 
vriQ = 1.5 TeV (signal), superimposed on the expected Standard Model background (SM), 
for 100 fb _1 of integrated luminosity. The mass window used to select the signal is indicated 
by arrows. 



We now consider the model-independent detection limit for such a resonance. 
The minimum number of signal events, iV™ n , needed to detect the resonance above 
the background is taken to be 5\/N b or 10 events, whichever is greater, in order to 
ensure that the signal is statistically significant. The minimum cross-section times 
branching ratio required to produce this number of events is also given in Table [3], 
after correcting for the detector acceptance and efficiency, and assuming an integrated 
luminosity of 100 fb -1 . 

It can be seen from these results that the background becomes negligible for 
graviton masses above 1800 GeV. It should be emphasised that these results are 
independent of the test model, and rely only on the assumption that the width of 
the resonance to be detected is negligible compared to the detector mass resolution. 
However, as an example, we can compare this cross-section limit with the values 
predicted by our test model, in order to determine the physics reach for ATLAS for 
this specific case. This is done in Figure |3], which displays a ■ B for the test model 
and (a ■ B) mtn . The test model could be detected for graviton masses up to 2080 
GeV. Using Figure |3|, similar estimates of the ATLAS search reach for other models 
can be determined without repeating the detector simulation or background analysis. 
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The search reach for other collider experiments, such as CMS, would require a new 
simulation with the appropriate mass resolution. 

Our results differ from those in [|b| for the case of a Z' with Standard Model 
couplings, both because of the acceptance difference arising from the spin of the 
resonance, and the Z' width, which exceeds 50 GeV for a Z' mass of 2 TeV. 
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▲ Test model cross section 
• Cross section limit 




500 1000 1500 2000 

Graviton mass (GeV) 

Figure 3: The cross-section times branching ratio, a ■ B, for G — ► e + e~ in the test model 
and the smallest detectable cross-section times branching ratio, (<r • B) min 



4. Angular distributions 

In order to demonstrate that the observed resonance is a graviton and not due to a 
spin-1 Z' , or similar exotic object, it is necessary to show that it is produced by a 
spin-2 intermediate state. The strong differences between the angular distributions 
shown in Figure [I] would allow the spin of the resonance to be determined given 
sufficient statistics. Information can also be gained on the relative production rate 
from the gg and qq processes. 

Figure f| shows the angular distribution expected in the test model for a graviton 
resonance mass of 1.5 TeV, after passing the events through the detector simulation 
and analysis cuts, using the statistics for 100 fb" 1 of data. The contributions from 
Standard Model Background, gg production and qq production of the graviton reso- 
nance are shown. In this case, 83% of the graviton production is from the gg process. 
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The distribution expected from a spin-1 resonance is also shown. The cutoff in the 
detector acceptance at \r)\ = 2.5 removes events at large |cos#*|. For heavy gravi- 
tons, which are produced with little longitudinal momentum, the effect is relatively 
sharp in cos#*, while for lighter gravitons and Drell-Yan processes, the acceptance 
loss reaches to lower I cos 9* I values. 




0.5 

cos(0*) 



Figure 4: The angular distribution of data (points with errors) in the test model for 
mo = 1.5 TeV and 100 fb _1 of integrated luminosity. The stacked histograms show the 
contributions from the Standard Model (SM), gg production (gg) and qq production (qq). 
The curve shows the distribution expected from a spin-1 resonance. 

A likelihood function was constructed to quantify the information in the angular 
distributions, defined as 

L = x q ■ f q (9*) ■ A q (M, 9*)/I q (M) + x g ■ f g (9*) ■ A g (M, 9*)/I g (M) 

+x D y ■ f DY (9*) ■ A DY (M, 9*)/I DY (M) (4.1) 

where Xi is the fraction of the events from each contributing process, fi{9*) is the 
angular distribution of the process, Aj(M, 9*) is the acceptance of the detector as a 
function of the mass of the electron pair and 9*, and 

h{M) = fi(9*) ■ Ai(M, 9*) dcos9* (4.2) 

i = q,g,DY for the processes qq — > G, gg — ► G, and qq — > Z/^* respectively. Only 
the shape of the distribution is used in the statistical tests, and the coefficients x are 
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constrained such that 

x q + x g + x DY = 1 (4.3) 

In order to evaluate the discovery reach of the experiment, in terms of its ability 
to reveal the spin-2 nature of the resonance, the following procedure was followed, 
intended to mimic an ensemble of possible experimental runs: 

1. N events were generated, containing the expected number Npy background 
events, with a Poisson distribution, and the remainder being graviton resonance 
decays with the qq and gg mixture of the test model at that mass point. 

2. — logL was evaluated for two hypotheses: 

= 0, xdy = 1, corresponding to a spin-1 resonance, such as a Z' 

(b) all the coefficients set to the values in the test model with the expected 
Standard Model background fraction. 

3. N was increased by adding graviton resonance events, until the difference in 
— logL between the two hypotheses reached 1.34, then 1.92 and finally 3.32. 
At these points the spin-1 hypothesis can be eliminated with 90, 95 and 99% 
confidence respectively. 

4. This procedure was repeated for many sets of events, and the average value of 
N needed for discrimination was found for each confidence level. 

This value of N provides a measure of the statistics needed in a typical experiment 
in order to determine the spin-2 nature of the resonance. N™ in = N — Npy is the 
number of signal events in the sample at the confidence limit. Table |4] shows results 
of this procedure. 

The 90% confidence limit for me = 500 GeV is raised by the requirement that 
the sample contains at least h^/Wjyy signal events. In this case, less events would be 
sufficient to establish the spin of the resonance, but would not meet our criterion for 
a discovery. At very high ma, the discrimination becomes more difficult because of 
the increased fraction of the resonance production coming from the qq initial state. 

Figure [5] shows the confidence limit curves and the corresponding cross-section 
times branching ratio in our test model. Smooth curves have been drawn through 
the points since some statistical fluctuations remain at high mass. These results 
are model independent, as long as the assumption that the coupling of the graviton 
resonance is univeral holds true. This assumption affects the fraction of graviton 
resonance production from gg and qq initial states, and hence the shape of the 
angular distribution. The cross-section times branching ratio of the test model is 
shown for comparison. In this model, the spin of the resonance could be determined 
with 90% confidence for graviton resonance masses as high as 1720 GeV. 
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m G 




N min 






a ■ B (fb) 




(GeV) 


90% CL 


95% CL 


99% CL 


90% CL 


95% CL 


99% CL 


500 


140 


141 


226 


1.90 


1.92 


3.07 


1000 


59 


70 


99 


0.80 


0.94 


1.33 


1500 


43 


48 


70 


0.61 


0.68 


1.00 


1700 


41 


48 


65 


0.61 


0.71 


0.96 


1800 


29 


33 


58 


0.43 


0.49 


0.87 


1900 


32 


40 


64 


0.49 


0.61 


0.98 


2000 


36 


41 


69 


0.56 


0.64 


1.08 


2100 


31 


45 


59 


0.49 


0.71 


0.93 


2200 


29 


33 


55 


0.45 


0.52 


0.86 



Table 4: The minimum numbers of signal events, iV™ m , needed to distinguish the spin-2 
and spin-1 hypotheses at 90, 95 and 99% confidence, and the corresponding cross-section 
times branching ratio, for a run of 100 fb -1 . 

1 
PQ 
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10 2 
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-1 

10 

500 1000 1500 2000 

Graviton mass (GeV) 

Figure 5: The smallest cross-section times branching ratio for which the spin-2 hypothesis 
is favoured over the spin-1 hypothesis at 90, 95 and 99% confidence. The cross-section times 
branching ratio of the test model is also shown. 

5. Conclusions 

The discovery limits for the detection of the decay mode G — > e + e _ , derived above, 
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▲ Test model cross section 




show that the resonance can be detected up to a graviton resonance mass of 2080 
GeV in our test model. The limits are model independent as long as the graviton 
couplings are universal and give rise to narrow resonances, with widths less than the 
experimental resolution. The angular distribution of the lepton pair can be used to 
determine the spin of the state. In our test model, the angular distribution favours a 
spin-2 hypothesis over a spin-1 hypothesis at 90% confidence for graviton masses up 
to 1720 GeV. Once the spin of the resonance has been determined, it will be useful to 
check the universality of its couplings Hl4j| . Muon and jet production cross-sections 
in the resonance region as well as decays into photons, massive vector bosons, higgs 
bosons and top quarks could be used to check this universality. The bosonic final 
states would provide a further test of the spin, since their angular distributions will 
be different from that of the electrons considered above. They would also serve to 
rule out other models of high mass resonances. For example, the Z' will not couple 
to a two photon final state. 
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